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cn ■ Abstract 
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At the LHC, qLQL production is one of the main SUSY production processes in the MSSM, which 
occurs due to the chirality flip caused by the Majorana gluino mass. This process is one of the 
sources of same sign two lepton (SS2Z) events, however, gluino production also contributes to this 
channel. In this paper, we develop a method to identify gluino and squark production separately 
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Qh! in the SS2Z channel, based on cuts on the kinematical configuration of the jets. We applied the 

method to the MSSM, a model with an extended gluino sector, and the Littlest Higgs model with 



T-parity (LHT), and found a distinctive difference between these models when considering the 



^ . numbers of SS2Z events selected by the cuts. 
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I. INTRODUCTION 



The Standard Model (SM) describes interactions among elementary particles very well. 
However astrophysical and cosmological observations such as WMAP have confirmed the 
existence of dark matter (DM) that cannot be explained in the SM [1]. The particle which 
DM consists of does not leave detectable signals in the detectors at high energy collider 
experiments, because it should be weakly interacting. 

If new particles are produced at collider experiments and decay into visible particles 
and a DM particle, we can observe large missing transverse momentum (J$t) in the events. 
Many models which predict DM candidates have been proposed. Among them, the mini- 
mal supersymmetric standard model with conserved R-parity (MSSM) is an attractive one. 
The lightest supersymmetric particle (LSP) is stable and a DM candidate. In the MSSM, 
quadratic divergences in the Higgs mass radiative corrections cancel each other, therefore 
the fine tuning problem is solved. At the large hadron collider (LHC), a pp collider with 
^fs = 14 TeV starting its operation in 2007 0,1^] , discovery of the squarks q and gluino g - 
super partners of quarks and gluon - is possible for the masses up to 2.5 TeV by looking for 
an excess of events with large f$T- We can also measure their mass spectrum by studying 
the decay kinematics of the quarks and gluino decay chains if there are enough events. 

However, the discovery of a I$t signature does not necessarily mean the confirmation 
of the existence of supersymmetry. A similar mass spectrum and decay pattern might be 
obtained for the Universal Extra Dimension model (UED) \m and the Littlest Higgs model 

n n q 

with T parity (LHT) |5|, |6|, |7fl. To study the origin of the If? signature, therefore, it is 
important to measure the other features that are characteristic of the MSSM. 

Many analyses have already been carried out in this direction. Recentfvprocesses sensi- 
tive to the gqq Yukawa type coupling constant have been investigated in [8). The Yukawa 
type coupling constant for the gqq vertex is the same as the gauge coupling constant due 
to SUSY. They study the same sign two isolated-lepton (SS2Z) channel to estimate the pro- 
duction cross section oi^q^qi). The SS2Z channel is one of the major discovery channels 
for supersymmetry, which is studied in jt], Q, 11]. The process is in principle sensitive to 



the coupling, however they found a large background from gq production. To measure the 
coupling constant, it is important to measure the cross sections of sparticle production pro- 
cesses separately. In this paper, we also focus on the cr(q L q L ). We will give a new method 
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based on cuts on the numbers of jets in the hemispheres for the purpose of separating q^L 
production from gq~i production. 

The production process occurs through a chirality flip caused by the gluino majorana 
mass term m g gg. To study the sensitivity to the majorana nature of the gluino mass, we 
consider a model with an extended gluino sector. This extension is inspired by the model 
that extends SUSY to iV = 2 in 12, In this model, an adjoint matter a is 

introduced, then the gluino can have a Dirac mass term m^ga. The gluino mass receives a 
contribution from mp. We also discuss the Littlest Higgs model with T-parity (LHT). This 
model contains quark partners and gauge boson partners (Wh, Ah) which decay into 
the stable lightest T-odd particle (LTP) Ah and SM particles as in the MSSM. Then, collider 
signatures are similar. There is progress concerning spin studies at LHC to distinguish these 
models flf] . \v\ . In this paper, however, we focus only on the difference among the production 
cross sections in these models. 

This paper is organized as follows. In Section 2, we first discuss the mass dependencies of 
the production cross section of q and g in the MSSM. The ulUl production cross section is 
typically of the order of 100 fb at trig ~ m g ~ 1 TeV, which may be detectable at LHC. The 
production processes of q^L and quqn occur due to the gluino majorana mass, therefore the 
mass dependencies are different from those of cx(<?L<?ij) and &(qLq* L ). We compare ct^^l) 
with a(gqi) and find that cr(q L q L ) <C <j(Ml)- In this section, we also choose a few model 
points for later analyses. 

For the MSSM with an extended gluino sector, the g, q production cross sections are 
functions of the three mass parameters - two majorana masses m g , itla and one Dirac mass 
ttid, where the mass terms are of the form m g gg + m^aa + mo{ag + ga). In particular, 
<y{qLq~L) becomes zero for some particular choice of parameters. We also discuss the LHT 
model. The production process occursby heavy SU(2) gauge boson exchange. <r(g_g_) 
is large compared with a{qq) of the MSSM [18]. On the other hand, because no gluon partner 
exists in this model, there is no problematic background corresponding to the gq production 
of the MSSM. 

In Section 3, We study SS2Z events to estimate cr(q L q L ) in the MSSM. To reduce SS2/ 
events from gqi, we use a fc -jet veto and the numbers of jets in hemispheres defined following 
the procedure proposed in [19j. We demonstrate that production processes gg, gq, qq can be 
distinguished by the efficiency under the cuts. We also discuss the dominant ti background. 
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In Section 4, we calculate the expected number of SS2Z events in the model with an 
extended gluino sector considered in Section 2. We discuss the sensitivity of the QlQl pro- 
duction cross section to the majorana gluino mass m g at a few model points. We also 
estimate the number of SS2Z signature in the LHT model. The efficiency under the cuts 
to reduce gluino background turns out to be useful to prove the existence of quark partner 
productions and difference from the MSSM prediction. Section 5 is devoted to the summary. 



II. PRODUCTION CROSS SECTIONS AT LHC 
A. MSSM production cross sections at LHC 

In the MSSM, sparticles are always pair produced at LHC because of R-parity conserva- 
tion. Production processes gg, gq, qq occur copiously unless the masses are much heavier 
than 1 TeV. Because u quark and d quark parton distribution functions (PDF) of a pro- 
ton are much harder than the other partons, u and d are mainly produced among squarks. 
In particular, production processes ulUl, d L dL, urUr, dudji, etc. require a chirality flip, 
therefore they do not occur if the gluino mass is not of majorana type m g \\. Only gluino 
exchange diagrams contribute to the productions. 1 

If the sparticle mass spectrum is known, the production cross sections at LHC are cal- 
culable. The mass dependencies of some sparticle production cross sections are shown in 



Figures HH3j Here, we use CTEQ 61 20| as the PDF. The horizontal axis is the gluino 
mass m g , and the vertical axis is the squark mass rriq. 2 The production cross section of each 
process is shown in contour lines in units of pb. 

^(Q'-L^l) (<?l = ul or di) are shown in Figures [1] a ~ c. <t(ulUl) is 0.05 pb at rriq = m g = 
1000 GeV, and 0.12 pb at rriq = m g = 800 GeV. This shows that they strongly depend on the 
squark mass rriq. The mass measurement error at LHC for squarks and gluino can be around 
a few percent, if the number of events is large enough 21|. u{qLqL) changes by 10% when rriq 



changes by 3% around rriq = m g = 1000 GeV and rriq = m g = 800 GeV. On the other hand, 
the cross section only weakly depends on gluino mass. This is because the amplitude has a 



1 Here, we neglect the contributions from neutralino and chargino exchange diagrams. This assumption is 
reasonable, because g 2 < g 2 ,m^r < nig in mSUGRA. 

2 We set niu L = m q , mg = rriq + 6 GeV, ra iLR = rriq — 19 GeV. 
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FIG. 1: Contour plot of a) (j{ulul)-, b) cr^d^di), c) aiu^di) and d) o{ulur) as a function of m g 
and ra„. 
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FIG. 2: Contour plot of a) a{gui) and b) o{gdjJ) as a function of m g and m ? -. 



factor of rrig from the chirality flip, which compensates the suppression from the t-channel 
propagator. <j(ulUl) changes by up to 10% in the range 500 GeV< m g < 1100 GeV around 
rriq = 800 GeV, and also changes by up to 10% in the range 700 GeV< m g < 1300 GeV 
around rriq = 1000 GeV (Figure [Ih.). In the case of m g = 3TeV, a{uL,Ui) is 0.04 pb at 



800 GeV, and 0.02 pb at rriq = 1000 GeV. As a result, even when the accuracy of 



the mass measurement of the gluino is bad, the qi^L production cross section is a useful 



5 



m q (GeV) 



0001 








- 


0.003 








- 


-0,005^ 
















* D 




0.03 
0.05 "' 


" A 


9 

















600 800 1000 1200 1400 1600 1800 2000 
m g (GeV) 

FIG. 3: Contour plot of a(uLd* L ) as a function of m g and rriq. 

observable that can be used to quantitatively test the MSSM. The behavior of the qnqn 
production cross section is the same. 

^(QlQr) (such as ct^ulUji) etc.) depends on the gluino mass more sensitively than cr(q L q L ). 
It decreases as the gluino mass increases. <j{ulUr) drops by half as m g increases from 
500 GeV to 1100 GeV for rrig = 800 GeV, and also drops by half as m g increases from 
700 GeV to 1300 GeV for nig = 1000 GeV (Figure [Hi). cr(gq L ) depends on the gluino mass 
even more because an onshell gluino has to be produced (Figure |2^j2b). 

By investigating q^qh production processes, we can probe the majorana nature of gluino 
mass. To measure the production cross section of ulUl and d^L, same sign two lepton 
(SS2Z) events are useful, which is studied in [8]. This idea is as follows: a Wl(Jl) may 
dominantly produces / + (/~) through the decays of 

ul -> xtd^l+Vid -> x°il + vid d L -> Xi u -> l~ViU -> y§TviU 

or -^x\W + d -> xil + "id , or -> x\W~u -> jftTvw . (1) 

Therefore l + l + events are sensitive to UlUl production and events to dj^di production. 
This signature implies the existence of the Yukawa type vertex ggq. 

The ratio of l + l + and has more information, as studied in [llj]. The ratio of the 
fractions of u and d in the PDF is 2:1. Then ct(ulUl) is larger than c^g^o^) and their ratio 
is about 4:1 (Figure HK,b). Thus the ratio of the number of SS2Z events N(l + l + from UlUl) '■ 
N(l~l~ from d^di) should be 4:1 if leptonic branching ratios of ul and di are the same. 

The processes involving u* L and d* L etc. also become sources of SS2/ events although the 
cross section is not large (a(uLd* L ) ~ 0.1cr(w£,tti) in Figure [3]). Since the basic observable 
is the sign of the leptons, in the following we define q[ as {ul, d* L , cl, §* l } which can be a 
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parent of Z + , and ql as {u* L , cIl, c* l , Sl} which can be a parent of l~ . 

The mass dependencies of the production cross sections of q^, ql are shown in Figures 
HMHJ By comparing Figure H^i with Figure [1^, and Figure [4b with Figure [lb, we can see 
that UlUl (dLd L ) is dominant in q[q[ (qlql) respectively. By comparing Figured with 
Figure Eh, and Figure [5b with Figure Eb, we also find that qul (gdi), is dominant in gcfc 
(Ml) respectively. 

The g can decay into qfq T , therefore gq^ production also produces l ± l ± . Moreover, 
cr(gqf) is larger than cr(qfqf) unless the gluino is too heavy (Typically u^gq^) ~ 5a(g ± g ± ) 
), thus the gq^ production process becomes background to the q^q^ production process. 
The ratio cr(gq + ) '■ cr(gq~) is about 2:1 (Figure [5]). Then the ratio of the numbers of SS2/ 
events N(l + l + from gq + ) : N(l~~l~ from gq~) should be 2:1. 

The gg production process also produces l^l^ and N(l + l + from gg) : N(l~l~ from gg) 
should be 1:1. This process, however, does not produce problematic background because 
the SS2Z branching ratio of gg production is small, although a(gg) may be larger than a(qq) 
in the mSUGRA model. 
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FIG. 4: Contour plot of a) a(q^q^), b) a(q L q L ) and c) cr(q^q L ) as a function of rrig and nig. 
The ratio of a(qfqf) / cr(gq^) increases as the gluino gets heavier or the squark gets lighter 
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FIG. 5: Contour plot of a) cr(gq^) and b) a(gq L ) as a function of m g and rriq. 
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(Figure [6]). When m g ~ rriq the ratio is always less than 1/2 for < 1500 GeV. 



B. The MSSM model points 

The numbers of SS2Z events from gq^ production and from qfqf production depend on 
their decay patterns and the branching ratios. For the study in the following sections, we 
choose a few model points to fix the branching ratios. We consider the model points with 
m-g > rriq. If rriq > m g , both of q^ decay into g and this g can produce a lepton or an 
anti-lepton with the same probability, and we cannot distinguish q q and q + q~ . 

We take relatively heavy sparticle masses (rriq,m g > 800 GeV), because a(qfqf) is too 



small compared with a(gq^) if the masses are lighter. For example, at SPSla (mq,m g ~ 



600 GeV), which is a popular benchmark point defined in [22], cr(gq) is about ten times as 
large as cr(qq). We choose four points A~D shown in TableUso that cr(qtQt) / a (Mt) ~ 1/5, 
because we find that one can reduce gq background by factor of 1/10 by applying various 
cuts in Section 3. The points A~D and SPSla are also marked in Figure [THE We can see 
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in Figure [6] that the ratio cr(qt<li)/< J (Mi) * s about 1/6, 1/4, 1/3, 1/2, and 1/10 for Points 
A~D and SPSla respectively. 

Points A and B are mSUGRA points where (m ,mi) are (100 GeV,340 GeV) and 
(100 GeV, 450 GeV), A = 0, tan/5 = 10, sign \i > respectively. At these points, the 
mass difference of fx and Xi is small. This feature is favored to reduce the LSP abun- 
dance. Points C and D are the mSUGRA points where (mo, mi) are (370 GeV, 340 GeV) 
and (400 GeV, 450 GeV), A = 0,tan/5 = 10, sign fi > respectively, except the low energy 
gluino masses are heavier than the mSUGRA predictions by 300 GeV. By increasing the 
gluino mass, the production cross section of the gluino decreases, therefore gq backgrounds 
are smaller. Sleptons are heavier than the lighter charginos xt > thus the xf does not decay 
into a slepton but decays into a boson. These points are not favored cosmologically, 
however the mass density can be reduced by tuning the mass of the pseudoscalar higgs boson 
mp as mp ~ 2rriLSP without changing the rate of SS2/ signal if we go beyond mSUGRA. 
Moreover, the decay patterns have some similarity to those predicted in the LHT as we will 
see later. SPSla is also written in Table HI and HT1 for reference. 

The masses of some particles at our model points are shown in Table HI These spectra 



are calculated using ISA JET 23|, |24j. The other mass spectra are given in the Appendix 
IA 11 The SUSY production cross sections of our model points are also shown in Table [Til 



This is calculated with HERWIG 6.5 [25J, where we use the CTEQ 61 PDF 



20]. 





mo mi Aq 
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m-g mu L mu R m i± m~ hx mf m^o 


Point A 
Point B 
Point C 
Point D 
SPSla 


100 340 
100 450 
370 340 
400 450 
100 250 -100 


809.86 737.25 714.56 559.18 683.97 160.79 256.36 132.74 
1047.83 951.16 919.51 734.57 883.16 194.28 324.47 179.11 
1123.23 808.67 787.79 585.39 731.10 387.18 437.13 133.95 
1360.22 1017.91 988.10 804.20 928.37 429.40 503.38 180.54 

595.19 537.04 520.45 379.14 491.92 133.39 202.12 96.05 



TABLE I: The mass spectra at Point A~D and SPSla. Here the unit of masses is GeV. 



C. The model with an extended gluino sector 



In this section, we consider a model with an extended gluino sector. This model was 



originally considered in 
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15| as a model with enhanced particle contents of N = 2 





a(SUSY) 




a (Qr Qt ) 


a (qi<ir ) 




a(gq T ) a(gg) 


Point A 


8.621 


0.2251 


0.0674 


0.4247 


1.3580 


0.6005 0.7134 


Point B 


2.023 


0.0750 


0.0189 


0.1309 


0.2949 


0.1208 0.1385 


Point C 


3.418 


0.1321 


0.0371 


0.2572 


0.3410 


0.1406 0.0989 


Point D 


0.963 


0.0494 


0.0118 


0.0875 


0.0897 


0.0345 0.0280 


SPSla 


45.890 


0.8033 


0.2868 


1.600 


7.408 


3.544 4.872 



TABLE II: The SUSY production cross sections for some processes at Point A~D and SPSla. Here 
the unit of the cross sections is pb. 

SUSY to solve the little hierarchy problem. The model has new fermions af in adjoint 
representations of each SM gauge group Gj. The majorana gaugino A° has a Dirac mass 
term with a". 

Inspired by the model, we consider phenomenologies of the gluino sector. The mass term 
of the gluino is extended as follows, 



/'mass 

L, 9 



-m g gg 



1 f- _> 
9 5 




(2) 



Here, g and a are four component spinors that satisfy the majorana condition. The m g 
is a majorana mass for g, 3 mo is a Dirac mass between g and a, tua is a majorana mass for 
a. We leave the SU(2) and U(l) gaugino sectors unchanged. 

The mass eigenstates gi, c/2 are majorana particles, and the mass eigenvalues are given by 



m 



91,2 



i (^m g + m A ± yj (m g - m A ) 2 + 4m? D 



\ m 9i \ < \ m 92\- 



(3) 



The mass eigenstates are defined as follows, 




cos <p sin 
— sin 6 cos 




(tan0 



m gi - m g 
m D 



)• 



(4) 



\92 / 

In the limit of nip 0,rriA — > oo, a decouples from the MSSM fields, then the phe- 
nomenology becomes identical to that of the MSSM. In the limit of m g = 0, m A = 0,mD / 
(we call this limit the pure Dirac limit), the masses of the two gluinos become the same, 



3 In Rcf. 



13 ] . rrig is taken as zero. 
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then gi and cj2 interfere strongly so that g and a form a Dirac particle go and its anti par- 
ticle go- The model has continuous R-symmetry in the pure Dirac limit. We can assign 
R-charge 1 for g D and u L ,—1 for u R and so on. The pair production of UlUl, UrUr through 
gluino exchange, which is one of main production process in the MSSM at the LHC, is 
forbidden by the R-charge conservation law. This can be also understood by the fact that 
existence of nonzero majorana mass of gluino is necessary for q^L, QrQr processes. While 
QdUl pah production is allowed, (JdUl pair production is forbidden in the limit, and once 
go is produced, go can decay into ulu and did but can not decay into u* L u nor d* L d. The 
difference between the pure Dirac limit and the MSSM is clear. If ttia, m g <C m^, the mass 
difference of two gluinos Am^ 2 is small, continuous R-symmetry exists approximately, and 
the phenomenology is similar to that in the pure Dirac limit. We do not investigate the 
phenomenology of this case any further. 

This model has two gluino like particles. We may be able to observe two gluinos if the mass 
difference is large enough compared with the decay widths of the gluinos (r§ 12 <C Amj, 2 ). 
In this case, each gluino decay produces q + and q~ with the same branching ratio. If g<i is 
too heavy so that the production cross section is too small to be observed at LHC, we can 
observe only g\. It is not possible in this case to distinguish this model from the MSSM only 
by the mass spectrum. We focus on this case and study the deviation of the production 
cross sections from the MSSM predictions. The gluino sector of this model has two degrees 
of freedom in addition to m gi and we take them as the other gluino like particle's mass m g2 
and the majorana mass m g . 

In Figure [3, we show the production cross sections of SUSY processes such as q^q^ and 
gqf as functions of the majorana gluino mass m g . Here, we fix m g2 = —3000 GeV and 
m gi = 1047.83 GeV, which is the gluino mass at Point B. The c/2 cannot be searched for at 
LHC. We take the mass spectrum of the other sparticles to be the same as that of Point B. 

In the limit of m g = m gi = 1047.83 GeV, mo is zero and m g2 = rriA = —3000 GeV. The 
c/2 decouples from q, and the cross sections involving gi are the same as those of Point B. 
Changing m g from this value distorts the model from the MSSM. As the majorana mass m g 
decreases keeping m gi and m g2 fixed, the total SUSY production cross section decreases. In 
particular, cr(g^g^) decreases quickly, while ^{giq^) decreases linearly because cx(<7i , ? ± ) is 
proportional to cos 2 </> and m g = cos 2 4>(m gi — m g2 ) + m g2 . The fraction of g in g2 increases, 
but c/2 is too heavy so that cr(^2?) is small. a(gigi) is less sensitive to the majorana mass 
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m g . This means that the squark exchange diagram does not contribute much to the cross 
section. 

When the majorana mass m g = 0, cr(g^g^) is still nonzero. This is because the majorana 
mass of the adjoint fermion rriA causes the chirality flip. There are also minor contributions 
from a(uLd* L ) and so on, which are not suppressed by the m g factor. 

As m g decreases further (absolute value \m g \ increases), a^qj^qjj) approaches zero around 
m g ~ —1500 GeV, while cr(giq^) is reduced by factor of 3. This behavior can be explained 
by the dependence of the subprocess cross section a(qq — > qh^h) on the mass parameters as 
follows, 

& f 1 m m 4 ^' 



Pf f 

a(qq -> Ml) = ^ J d(cos6) 



oc 



9 

X \ T ij T kl{f( m 9l> m q: P) COs2 + f( m 92, m q: P) SU1 2 0} 

+TSTkj{f( m gi> ™» -P) C0S V + f( m §2, m q, -P) sin 2 0}| 
m gi cos 2 m g2 sin 2 . 



(5) 



where, f(mi,mf,p) = — = = — — ^ — — - — ; ~ — = — - — (6) 

v r ' m 2 + m 2 + 2|p| 2 -2£|p| cos 9 m? + mj w 

Here a(qq — > qhqh) is the cross section of the subprocess qq — > qiqL with center of mass 
energy s. p is the momentum of one of the created q^. We set the z-axis along the momentum 
of one of the initial quarks, 9 is the polar angle of p from the z-axis, and f3f is the beta factor 



(Pf — Jl — 4m 2 - 1 s). T a denotes a generator of the SU(3) group and g s denotes the gauge 
coupling. is the mixing angle between the g and a as defined in In the last line in 
eq.flH]) and the last equality in eq.flEJ), we take the limit that |p| <C rriq because a(pp — > qhqi) 
is dominated by the threshold production. We can calculate that crfe^) approaches zero 
around (f> = </> defined as 



ro*. ml + m 2 



tan o = A /-^l * \ . (7) 
m g2 m z gi + m~ 



This corresponds to m g = —1516 GeV for Point B (here, we use the relation m g = cos 2 <f)m gi + 
sm 2 (j)m g2 ). Note that we again neglect the contributions from chargino and neutralino 
exchange diagrams in the calculation. They are at most of the order of 10~ 3 pb, and 
negligible. 

<j (Ql^l)/ '(Ml) i s l ess than 5% in the range of —1800 GeV < m g < —800 GeV, while 
^(^lQ'l) / /(X (Ml) * s aDou t 25% in the MSSM limit {m g = 1047.83 GeV). The production cross 

12 



section <j{qLQL) is reduced by more than factor of 5 compared with gqi in this range (Figure 



As m g decreases further, cr(q L q L ) increases again while a(giq L ) keeps decreasing. When 
= —3000 GeV, cj2 is g and g\ is a. In this limit, c(g L g L ) is enhanced by the factor of m g 
in the amplitude, and cr(giqf) = because g\ does not couple to q. We do not discuss this 
region because the production and decay pattern would be significantly different from the 
MSSM. 





-3000 -2500 -2000 



-1500 -1000 
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500 1000 



FIG. 7: The production cross sections as a 
function of the gluino majorana mass m g for 
the model with an extended gluino sector. The 
mass spectrum is the same as that at Point B. 



FIG. 8: The ratio of ^{q^q^) and &(Ml) as a 
function of the gluino majorana mass m g for 
the model with an extended gluino sector. The 
mass spectrum is the same as that at Point B. 



Note that the mass spectrum is the same as that at Point B for entire m g in Figure [71 We 
can detect the deviation from the MSSM only through measurements of production cross 
sections such as a(qfqf), a(gqf) or their ratios. 



D. The Littlest Higgs model with T-parity 

The Littlest Higgs model with T-parity (LHT) 

Baa is an alternative scenario that 

solves the quadratic divergence problem for Higgs mass andpredicts a stable DM candidate. 
This model is the extension of the Littlest Higgs model [26[]. The features of the Littlest 
Higgs model are: 

1. Higgs bosons are introduced as pseudo Nambu-Goldstone bosons of global symmetry 
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breaking SU(5)/SO(5). The global symmetry is partially gauged, and the gauge 
symmetry is [577(2) x U(l)] 2 . 

2. The gauge groups are spontaneously broken at scale / as [SU (2) x U{1)} 2 — > SU (2) L x 
U(1)y- Heavy gauge bosons of the broken gauge groups are called as W H , Zh, Ah- 
The masses are, 



M Zh ~ M Wh ~ gf 



1-^ 

8P 



g'f 



5v 



8/ 2 

Here, w is the electroweak symmetry breaking scale, v ~ 246 GeV. 



(8) 



3. To cancel the quadratic divergence of Higgs mass, the third generation fermion sector 
has to be extended to respect the global symmetry of the theory. In particular, T + 
has to be introduced as a partner of the t quark. The T + is a SU(2)l singlet Dirac 
fermion. 

However, this model suffers from large tree level corrections to the electroweak parameters. 



Even if the parameters of the models are tuned to reduce 



;he corrections, / becomes large 



enough that the fine tuning problem is reintroduced |27l . |28| . To solve this problem, the 
LHT model imposes invariance under the T-parity that corresponds to switching the two 
[SU(2)xU(l)] gauge groups. Matter sectors are extended so that there is a T-odd partner 
for each SM fermion. 

The T-parity plays a similar role to the R-parity of the SUSY model. The Lightest T- 
parity odd particle (LTP) cannot decay because T-parity is multiplicatively conserved for all 
vertices. The LTP can be a candidate for DM. Moreover, T-odd particles can be produced 
only in pairs in a collider experiment and each of them must decay into final states including 
an odd number of T-odd particles. As a result, the final states include at least two LTPs. 
The collider signal at LHC is large transverse missing energy I$t, just like SUSY. 

This model predicts a set of new particles. Amongst them, the heavy gauge bosons 
Wjji Zh, Ah, the SM fermion partners w_, gL, s_, c_, 6_, and the top partner T_ are T- 
odd. In the following, gl denotes {u-, c_, gL, s_, qZ denotes {gL, s_, 6_, U-, c_, £_}, 
g_ denotes qZ and qZ- 

No T-odd partner for the SU(3) gauge boson is introduced. On the other hand, the 



decay pattern of T-odd g_ is similar to that of q in SUSY. According to Ref. 



1811. about 



60% of q± decays into W^, 100% of Wg decays into W + and 25% of W + decays into 
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leptons. Therefore about 15% of q_ decay leptonically. This decay pattern is similar to 
that of qL at Point C. Thus, for the LHT model, we obtain SUSY-like signal as if there 
are only qq and qq* production at LHC Although there is no t-channel colored particle 
exchange, 0"(g_g_) is non zero due to the t-channel exchange of Wh and Ajj. It is higher 
than that of the MSSM because g_'s are fermions. For example, a{q J Zq J Z) is 0.7 pb, u{qZqZ) 
is 0.15 pb at M q _ = 800 GeV and / = 560 GeV, and cr(q±q±) is 0.2 pb, a(qZqZ) is 0.04 pb at 
M q _ = 1000 GeV and / = 700 GeV. <j{q±q±) is 4 ~ 5 times as large as a{q + q + ) at the MSSM 
model points with the same mass scale (See FigureHk). The ratio cr(qZqZ)/<j(qZqZ) is similar 
to a{q + q + )/a{q-q-) of the MSSM. The ratio a(l + l + )/a(l~l-) of the LHT is expected to be 
higher than that of the MSSM, because gq production dominates the total SUSY production. 
As we will see later, the production cross sections and their ratios will help to distinguish 
the LHT and the MSSM. 

E. Summary of the production cross sections 

There are several models which predict a MSSM-like collider signature with large If*?- 
In section 2.2 we have shown that the qiqi production cross section changes significantly if 
the majorana mass contribution to the gluino mass is reduced in the model with an extended 
gluino sector. In particular, in the case of m g < m§ u uiqiqi) can be reduced significantly 
compared to cr^gqi). 

On the other hand, for the LHT model, g_ may have a similar decay pattern to q while 
there is no particle corresponding to gluino. Hence, the signal is similar than that of the 
MSSM with an undetectably heavy gluino. <r(g_g_) is larger than cr(qq) by a factor of 5. 

The signal cross sections of these models are different to those of the MSSM, even if the 
mass spectrum is the same. Therefore u{qLqL —> l^l + X) is one of the key observables for 
MSSM studies. 

III. SEPARATION OF qq AND gq PRODUCTIONS 
A. Branching ratios of q, g 

To identify q^L production (mainly UlUl production), q^L —> 2/ + X events would be 
useful because BK(q£ -> Z± + X) > BK(qf ^ l T + X). In this paper, we assume 100% 



15 



of qn decays into q and the LSP Xv This is indeed realized over most of the mSUGRA 
parameter space. However, in less constrained models neutralino mixing can be different, 
such that also q~R has cascade decays involving the heavier neutralinos and charginos. In this 
case, q~L<iR and qnqn production can both produce SS2Z events, and the following analyses 
are more complicated. In particular, <j{qLqR) is not so sensitive to the majorana nature of 
gluino, therefore it is more difficult to probe the majorana nature of the gluino using the 
SS2Z channel. We assume that the LSP Xi is dominantly a Bino and the x® a Wino, so q~R 
does not produce leptons. 

The signal rate depends on the leptonic branching ratios of sparticles. The branching 
ratios at Points A~D are shown in Appendix IA2[ We summarize them in Table IIII1 
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TABLE III: Branching ratios of squarks 
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and gluinos at Point A~D. These are 
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We can see that qf produces when q^ decays through a chargino. BTZ(qf — > / ± ) is 
about 30% at Point A. If qf decays through a neutralino, q^ may also produce P. 

On the other hand, a gluino decays into third generation squarks q% more than the 1st 
and 2nd generation squarks g^, and also decays into qn more than qi at these points, 
because gluino and squark masses are close and the phase space of the gluino decay is 
sensitive to the small differences of squark masses. (The third generation squarks t, b are 
lighter than the 1st, 2nd generation q^, and q~R is lighter than q^.) We also find that 
B1Z(g — > l + ) = B1Z(g — > and they are small (8%), because g decays dominantly into <jr. 
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As a result, BTl{qtqt -> l + l + ) is 9%, i3ft(^ -> is 2.4%, Bft(<7£ -> 1+1+) is 0.6% at 

Point A. 

To measure cr(g^g^), we need to reduce the events involving g by means of appropriate 
cuts. We can achieve this in part by rejecting events with 6-tagged jets because BTZ(g — ► t 
or b) is large; BlZ(Jj — > I without 6-quark) is only 6%. For pair production processes, 
BTZ{qlql -> /+/+ without 6-quark) is ~ 9%, BTZ(gqt ~> without 6-quark) ~ 1%, 

B7Z(gg — > / + / + without 6-quark) ~ 0.1%. However, the efficiency of 6-veto is at most 60%. 

In the following, we neglect SS2Z events from gg production. This is reasonable if squark 
and gluino masses are sufficiently large, such as at Points A ~ D, because <j(gg) is small 
(Table HI]) and BTZ(g — > q£) <^ 1. Furthermore, the events contain more 6-jets on average. 



B. Event generation and detector simulation 



We generate about 300,000 SUSY events using HERWIG 6.5 [25j at Points A~D and 
SPSla. The number of events actually produced by HERWIG and the corresponding 
integrated luminosities are listed in Table HVl N(q^q^) : N(q£q£) is about 4:1 and 
N{Ml) • N(gq"[) is about 2:1 for these model points as discussed in Sec III Al 





N(SUSY) N{qtq + L ) 




N(qlq~ L ) N{M + L ) 


N{Ml) 


N(gg) 


jdtC 


Point A 


289906 


7865 


10698 


2197 


44007 


19330 


21065 


33.63 


Point B 


284544 


10601 


11115 


2526 


39208 


15818 


13242 


140.65 


Point C 


295042 


11411 


14666 


3072 


25793 


10729 


4748 


86.32 


Point D 


295695 


14505 


14394 


3279 


23589 


8982 


3333 


307.06 


SPSla 


293161 


5412 


10423 


1849 


46072 


22241 


31371 


6.39 



TABLE IV: Numbers of events generated by HERWIG 6.5. Here, the unit of the integrated lumi- 
nosities is fb _1 . 



We use AcerDET 29j for event reconstruction. AcerDET is a fast simulation and recon- 
struction package. It finds jets, isolated electrons, muons, photons and calculates the missing 
transverse energy from particles in the events. The granularity of the calorimetric cells is 
assumed as (0.1 x 0.1) in (rj x <fi) coordinates for \rj\ < 3.2. The clusters with pt > 15 GeV for 
Ai?conc = 0.4 are classified as jets. It also labels a jet as a 6-jet if a 6-quark with momentum 
Pt > 5 GeV is found within the cone AR = 0.2. The tagging efficiency of the algorithm is 
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about 80% and it is too high compared with the full simulation result of 60% in ATLAS. 
Therefore we assume that 60% of the 6-labeled jets are tagged. Isolation criteria for muons, 
electrons and photons are p? > 10 GeV and \rj\ < 2.5, separation by AR > 0.4 from other 
clusters and Ex < 10 GeV in a cone AR < 0.2 around them. For electrons and photons, 
we require Ai? cluster < 0.1. The E T is defined as 



*> T 



visible 



and calculated by summing the transverse momenta of all cells as follows, 



E 7 



cells 



(9) 



(10) 



The numbers of SS2Z events from each production process are shown in Table IVl for Point 



B. 



Point B 


generated 


1+1+ 


rr 


all cq 


Cl 


all 


C Cl 


total 


284544 


4363 2573 1680 


2231 1288 749 




10601 


1410 967 


958 


6 


2 2 




2526 


1 1 


1 


399 


264 258 


QlU 


11115 


88 54 


52 


112 


70 68 


Ml 


39208 


1720 1067 


467 


149 


84 31 


Ml 


15818 


46 31 


14 


732 


469 235 


M 


13242 


220 117 


22 


225 


121 26 



TABLE V: Numbers of SS2Z events from each production process at Point B 



Here, all, c and c\ denote the different cuts applied to the events, 
all : all SS2Z events generated by HERWIG 6.5. 

c : basic cuts, B T > 200 GeV, M cS > 500 GeV, B T > 0.2M eff and n 100 > 2 

(nioo is the number of jets with pt > 100 GeV) 
ci : Co and n& = 0. (n& is the number of 6-tagged jets.) 

M e ff is defined as 



(11) 



jets 
!p T |>50 GeV 



leptons 
|p T |>10GeV 
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The SS2/ events are mainly produced by q^q^ : gqti 99- 4 We find that N(l + l + from 
QlQl)'-N(1~1~ from q L q L ) is nearly 4:1 and N(l + l + from gq~^):N(l~l~ from gq^) is nearly 
2:1. Note that / ± / ± events are also produced from q + q~, if q decays into a neutralino and 
the neutralino decays into t + t~. After the Cq cut, contributions from qfqf,q^q2i Ml '99 
dominate the SS2/ events. We describe the number of generated l^l^ events after cut q as 
N^^; Cj). Comparing 7V(Z ± Z ± ; c ) with A(/ ± / ± ; ci), we see that the 6- veto cut (ci) reduces 
only events involving g. The l + l + events from gq^ are reduced by half. The Z + / + events 
from gg are reduced by one fifth. After the Ci cut, / + / + events from q + q + are 60% of all 
/ + / + events, while they are 40% under the c cut. 

C. Hemisphere cuts 

To further reduce events involving g, we next study the number of jets emitted from q 
and g. A gluino decays into a squark and a quark, and the squark decays into a chargino or 
a neutralino and a quark. Thus, a gluino usually emits at least two jets while a squark emits 
at least one jet. We can distinguish the parent particles by the number of high p T jets in the 
events. In this paper, we divide the final state particles into two groups called hemispheres, 
then we investigate the number of jets and the invariant masses in each hemisphere. 

SUSY production processes always occur in pairs due to R-parity conservation. Particles 
from each sparticle decay with momentum pi,pj,... are kinematically constrained so that 
(Pi + Pj + •••) 2 = m2 where m is the mass of the parent sparticle. When the parent sparticle 
is boosted the decay products are boosted in the same direction. We therefore divide all 
high p T objects into two groups: hemisphere 1 {pik}, hemisphere 2 {p2k}, which satisfy the 
following conditions. 

Any pu G {pik}, P2i £ {P2k} satisfy the conditions 
d(pi, ax ,pii) < d(p 2 , ax ,pii), 

d(p2, ax ,P2i) < d(p 1>ax ,p 2i ). (12) 

We define the axsis of hemisphere p ax and the distance between two 4-vectors d(pi,p2) 

4 The other SS2Z events come mainly from productions involving third generation squarks, charginos and 
gluino. 



19 



as follows. 



ax — 




(13) 



^(Pax, Pi) 



ax 



Pax | COS 9i)E ; 



ax 



(Here, 0j is the angle between p ax and pj). (14) 



(£ ax + £,) 2 



Here, high objects mean jets with > 50 GeV and r] < 3, leptons, photons with 
p T > 10 GeV and r] < 2.5. 

Our algorithm to find the hemisphere axes is as follows. We take the highest pt object 
in the event and regard its momentum as pi jax . Next, P2,ax is taken as the momentum 
of the object which has largest |p|Ai2, where p is the momentum of an object, AR = 
a/ (At]) 2 + (A0) 2 , and Arj and A0 are the differences of the pseudo-rapidity and azimuthal 
angle of p from pi )ax respectively. We calculate d{p\ i3X ,p) and d(p2,a, x ,p) for all high pt 
objects. We identify it to hemisphere 1 if d(pi )ax ,p) < d(p2,ax,p)- Otherwise, to hemisphere 
2. After that, we redefine pi t3x and P2,ax as the new hemisphere axes by using eq. f|T3|) . We 
iterate the same operation using new p l ax and p 2 , ax five times. 

After determination of two hemispheres, we calculate the maximum invariant mass rrijj 
of all jet pairs in a hemisphere. We call it rrijji for hemisphere 1 and rrijj2 for hemisphere 2. 
When a hemisphere has only one or zero jets, we define rrijj = 0. 

The 2-dim plots of rrijji vs. for production processes gg, gui, ulUl at Point B are 
shown in Figure [9k~c. Here, the plotted events are not only SS2Z events but all events after 
imposing c cut. 

We find gg events are mainly distributed in the region of rrijji 7^ and rrijji ^ 0, gq 
events are mainly distributed in the region of either rrijji = or rrijj2 = 0, qq events are 
mainly distributed in the region of rrijji = rrijj2 = 0. This is because g produces two high 
Pt jets while q produces only one high px jet. Therefore we require rrijji = m jj2 = as the 
cut to reduce g production events. 

D. Numerical results after hemisphere cuts 

We apply the cuts rrijji = and/or a b- veto on the SS2/ events and investigate the 
efficiencies for each production channel. We subsequently apply the cuts to remove the 



20 







a) gluino gluino 






1200 














70 




00 


D 70*1*04 


100 


1000 
300 






























to 

1 


600 














; 


400 
200 


: 










- 20 
















- 10 




200 400 600 800 1C 


00 1 200 



b) gluino u_L 



I.217E+05 31.0 



90C' 

80(' 

60C 
50C 

30C 

H20C' 
I0(i 



200 400 600 800 1 000 1 200 







c; 


u_L u_ 


L 


200 
















.00 




Z.00 


000 
300 




























600 
400 














200 














n 


■ 













2CO 



60C 
1V 

7'.' 



200 400 600 800 1000 1200 



FIG. 9: rrijjivs.mjj2 distribution for a) gg, b) guL, c) production events. 



contributions from gg production, while keeping the contributions from qq. We define cuts 
C2 ~ C5 in addition to Co and ci as follows, 

C2 : Co with min(mjji, m^) = 0. 

c 3 : c with m^i = 777,^2 = 0. 

C4 : Co, rib = and min(m--i, rrijj-i) = 0. 

C5 : Co, rib = and rrijji = rrijj2 = 0. 

ratio : A r (c 5 )/A r (c ). This describes the efficiency for c 5 cut. 

The number of the SS2Z events after these cuts are shown in the Table [VTl The cut requiring 
rrijji = or mjj 2 = (c 2 , c 4 ) reduces gg events drastically and gq events moderately. The 
cut requiring rrijji = and rnjj 2 = (c 3 , c 5 ) further reduces gq events. On the other hand, 
qq events survive under the cut C5 compared with gg, gq events. Note that a 6-veto cut is 
not essential to reduce the gluino contribution. The cut C3 reduces events involving gluino 
less than qiqL events although the 6-veto is not applied, gq events are dominant under the 
c cut and qq events are dominant under the c 5 cut. 

The ratio of N(l + l + ]C5)/N(l + l + ;co) for qq productions is more than 30% at our model 
points. It is 3 ~ 5% for gq production at Point A and B and about 1% at Point C and D. We 
can obtain a pure SS2Z event set from qq production by the C5 cut. Dominant contributions 
to N(l + l + from qq) are from qiqi production (mainly ulUl)- We can see in Appendix 
IA31 N(l + l + from u L u L and m l c l ) = 1256, N(l+l+ from u L d* L and u L s* L )=138, N(l+l+ from 
■Ui<ii)=84 are obtained among N(l + l + from all gg)=1479 at Point B. In Table |VH we also 
show the results at SPSla for a reference. The efficiencies of these cuts are similar to the 
other points. However contamination from gq production is larger. 
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Point A 


all 


Co 


Cl 


C2 


C3 


c 4 


C5 


ratio 


1+1+ total 


1519 


716 


349 


463 114 


267 


94 


0.131 


QQ 


235 


140 


132 


123 


57 


117 


56 


0.400 


9Q 


618 


333 


125 


205 


27 


83 


16 


0.048 


99 


184 


79 


14 


27 


2 


5 





0.000 


l~l~ total 


1213 


610 


286 


368 


89 


216 


75 


0.123 


QQ 


151 


92 


90 


83 


42 


81 


41 


0.446 


M 


472 


262 


108 


159 


19 


75 


13 


0.050 


99 


172 


93 


20 


31 


1 


6 





0.000 



Point B 


all 


co 


Cl 


C2 C 3 


c 4 c 5 


ratio 


1+1+ total 


4363 


2573 1612 1765 465 1293 441 


0.171 


QQ 


1479 


1014 1001 


894 363 


883 361 


0.356 


M 


1765 


1098 


433 


613 52 


279 36 


0.033 


99 


220 


117 


18 


26 1 


5 


0.000 


total 


2231 


1288 


708 


809 187 


561 171 


0.133 


QQ 


499 


326 


316 


286 122 


278 121 


0.371 


M 


861 


541 


237 


304 23 


169 15 


0.028 


99 


225 


121 


22 


32 


11 


0.000 



Experimentally, we can only observe the total number of SS2/ events. The ratio 
N(l + l + ; c 5 )/N(l + l + ; c ) for total events is 13~17% at our points. If there is no q L q L produc- 
tion, it becomes less than 5%. If the efficiency of N(l + l + ; c^,)/N(l + l + ; Cq) for qq and gq can 
be obtained by MC simulations and from the other constraints, we may estimate cr(qLqL) 
from N^fi). The model parameters we need are the gluino, squark masses to constrain the 
branching ratio of the gluino, leptonic branching ratios of squarks and their decay kinemat- 
ics. The ratio of events with 2 leptons to events with 1 lepton should be useful to estimate 
the leptonic branching ratio of q. The decay cascade g — > bb(tt) — > x tb(x + ^b) also emits 
leptons from W ± decays. The branching ratio must be estimated carefully from 6-tagged 
samples. Estimation of the errors on the branching ratio is beyond the scope of this paper. 
Information on the mass would be obtained from various end points of the decay distribu- 
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Point C 


all 


Co Ci c 2 c 3 c 4 c 5 


ratio 


l + l + total 


1081 


467 242 284 66 186 62 


0.133 


qq 


259 


154 144 133 57 127 57 


0.370 


M 


327 


165 55 69 3 29 2 


0.012 


99 


47 


10 1 2 


0.000 


total 


618 


233 110 134 20 83 17 


0.073 


qq 


77 


40 39 37 12 36 12 


0.300 


M 


164 


66 23 29 2 14 


0.000 


99 


32 


11 3 6 2 


0.000 



Point D 


all 


CO Cl c 2 c 3 c 4 c 5 


ratio 


l + l + total 


1157 


571 318 335 82 243 78 


0.137 


qq 


356 


216 207 182 72 179 72 


0.333 


M 


333 


169 56 64 2 26 2 


0.012 


99 


35 


18 5 7 4 


0.000 


l~l~ total 


588 


263 116 128 27 83 24 


0.091 


qq 


84 


55 47 45 16 40 16 


0.291 


M 


131 


63 25 23 1 14 1 


0.016 


99 


27 


9 


0.000 



spsla 


all 


CO Cl c 2 c 3 c 4 c 5 


ratio 


l + l + total 


1314 


461 250 326 80 193 63 


0.137 


qq 


164 


72 71 59 29 58 28 


0.389 


M 


550 


226 106 150 22 75 15 


0.066 


99 


155 


58 15 32 2 9 2 


0.034 


total 


1070 


345 180 242 64 147 46 


0.133 


qq 


109 


47 45 41 16 40 16 


0.340 


M 


394 


156 64 109 27 51 16 


0.103 


99 


147 


39 9 17 2 7 


0.000 



TABLE VI: Numbers of SS2Z events after the cuts for Point A~D and SPSla. Here, cq is set to 
$ T > 200 GeV, M cS > 500 GeV, tf T > 0.2M cS , n wo > 2. 
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tions 
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32 
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34 . |35| . In addition, there are various sources of systematic errors for 



the cross sections such as the PDF, NNLO corrections and so on. Following the arguments 
in j^l, we assume that the overall uncertainty on the total production cross section is 17%, 
which comes from 10% PDF uncertainty, 8% NNLO corrections [36] and 10% cr(?L?z) error 
from a 3% squark mass uncertainty. If the squark mass error is 10%, oici^qi) error is 40%, 
and this uncertainty dominates the systematic errors. 

Uncertainties on the cross sections may partly cancel in iV(Z ± / ± ; Cs)/A^(/ ± Z ± ; cq) ratios. 
The errors on the absolute sparticle masses are common for both a(q~L<lL) and &(gq~L), and 
rrig — uiq may be known precisely. If the cross section errors cancel, the total systematic 
error becomes 14%. Uncertainties on leptonic branching ratios of squarks may also partly 
cancel because the gluino emits leptons through its decay into a squark. 

Note that the ratio of the production cross sections of gui to gdi is 2:1, and that of it lUl 
to d^di is 4:1. Indeed, at Points C and D, the ratio N(l + l + ; C5) to N(l~l~; C5) is 4:1, while 
that is 2:1 for c . At Point A (B), the ratio of N(l + l + ; c 5 ) to N(l~l~; c 5 ) is 5:4 (5:2), while the 
ratio of iV(/ + / + ;co) to N(l~l~;co) is 7:6 (2:1). These discrepancies from &(q + q + )/a(q~q~) 
or a(gq + ) /cr{gq) m Table [IV] are caused by decays of the second lightest neutralinos from 
gj. According to Appendix IA 21 the left handed squark can decay into Xi T+T ~ through xt 
with a large branching ratio at Point A (B). Therefore there is a contamination from q£q^ 
(mainly Uhdi) production to SS2Z events and the ratio of N(l + l + ) to N(l~l~) gets closer to 
1:1. Therefore the charge of the hard lepton does not reflect the sign of the parent particle. 



E. Tight cut for heavy mass spectrum 

We now consider the SM background to SS2Z events. Here we only consider ti production 
which is found as dominant background in [8j . The background comes from the events where 
one lepton comes from the leptonic decay of a top quark while the other lepton comes from 
accidental sources such as b quark decays. The BTZ(t — > l~) is small, however the total 
number of ti events is significantly larger than the signal (cx(tt) = 400 pb at tree level). The 
ti events corresponding to 100fb _1 are also generated by HERWIG6.5, and the result is 
shown in Table IVIII 

At Points A and B, N(l+l+ from ti; c ) is less than 8% of N(l+l+ from SUSY;c ) and 
N(l+l+ from ti;c 5 ) is about 5% of N(l+l+ from SUSY;c 5 ). At Point C, it is 30% after c , 
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tt 


all 


C Ci c 2 c 3 c 4 c 5 


ratio 


l+l + 


1710 


152 70 133 32 65 17 


0.112 


rr 


1635 


146 76 133 20 70 10 


0.068 



TABLE VII: Numbers of SS2Z events from tt production after the cuts cq ~ C5 for JdtC = lOOfb . 
Here, cq ~ C5 are the same as Table 6. 

and 25% after C5. At Point D, they are ~ 1. In that case, more strict Iflx and M e g cuts are 
needed. 

In Table EH we have taken rf T > 200 GeV, M cff > 500 GeV and tf T > 0.2M cff as the 



Co cut. Numbers of ti after various I$t and M e g cuts are shown in Table IVIIII . The tt 
background is dramatically reduced if high Ifl T cuts are applied for large M e g [37(. Here, we 
change the c cut as follows. 
Ifi T and iVfeff cuts 



.(2) 



.(3) 



$ T > 200 GeV, tf T > 0.2M efr , M eS > 500 GeV 
$ T > 250 GeV, $ T > 0.2M efr , M cS > 750 GeV 
tf T > 300 GeV, $ T > 0.2M efr , M eS > 1000 GeV 



Number of high Pt jets 



ft-ioo > 2 : at least two jets with pt > 100 GeV 
^200 > 2 : at least two jets with p T > 200 GeV 



number of jets 


all 


^100 - 2 


™200 - 2 


tf T ,M eS cut 


all 


c (l) c (2) c (3) 


c (l) c (2) c (3) 


N(l+l+ from it; c ) 
N(l+l+ from ti;c 5 ) 


1710 

* 


152 38 9 
17 1 


13 6 4 
3 


N(l~l~ from ti;c ) 
N(l~l~ from tt;c 5 ) 


1635 

* 


146 43 15 
10 2 1 


11 8 2 




TABLE VIII: Number of SS2Z events from tt events after various Co cuts for 100 fb . 

We also show N^l^; c 5 )/A r (/ ± / ± ; c ) after the cut and the cuts on the number of 
high pt jets at Point D in Table HXl The numbers of events correspond to 307.06 pb^ 1 . We 
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can see that the efficiency iV(/ ± / ± ; C5)/A r (/ ± Z =t ; Co) for qq production and the efficiency for 
gq production only weakly depend on the basic cut and the cuts on the number of high 
Pt jets. The c 5 cut is still useful to reduce gq productions. 

We can drop the ti background without reducing the signal from SUSY events so much 
by taking and n 2 oo > 2 as Cq. However, it has recently been pointed out that the number 
of high pt jets increases significantly if matrix element (ME) corrections are included 38]. 
On the other hand, the I$t cut is not affected by ME corrections. The cut reduces 
background efficiently. 



number of jets 


all 


nioo > 2 


^200 > 2 


$T,M ef f cut 


all 


C (D 


c (2) 


c (3) 


c (l) 


c (2) c (3) 




N(l + l + from qq;co) 


356 


216 


211 


197 


162 


159 150 




N(l + l + from gq; Co) 


333 


169 


169 


161 


122 


122 119 




N(l + l + from qq;c§) 


* 


72 


67 


61 


59 


56 52 




N(l + l + from gq; C5) 


* 


2 


2 


2 


2 


2 2 


N(r 


H+;c 5 )/N(l+l+;co) for qq 


* 


0.333 0.318 0.310 


0.364 0.352 0.347 


N(l- 


-l + ;c 5 )/N(l + l+;c ) for gq 


* 


0.012 0.012 0.012 


0.016 0.016 0.017 




N(l~l~ from qq;co) 


84 


55 


54 


50 


38 


37 36 




N(l~l~ from gq; cq) 


131 


63 


60 


56 


47 


47 46 




N(l~l~ from qq;c§) 


* 


16 


15 


15 


14 


13 13 




N{l~l~ from gq; C5) 


* 


1 


1 


1 


1 


1 1 


N(r 


~r;c 5 )/iV(rZ - ;co) for qq 


* 


0.291 0.278 0.300 


0.368 0.351 0.361 


N(r 


~l~;c 5 )/N(l~l~',co) for gq 


* 


0.016 0.017 0.018 


0.021 0.021 0.022 



TABLE IX: Number of SS2Z events at Point D after various cq cuts for 307.06 fb 1 



F. Summary of the cuts 

It is important to reduce the background from gq production to measure the qLQL pro- 
duction cross section using SS2Z events. We give a systematic procedure to separate g and q 
based on the number of jets in a hemisphere, and demonstrate that it works well to separate 
qL<lL from gq for our model points. The hemisphere cut should work provided that rrig — nig 
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is sufficiently large that a jet from the decay g — > qq is detectable. 

Using the cut on the number of jets in a hemisphere and the fe-jet veto, SS2Z events from 
gq production are reduced by more than 95% while SS2Z events from qiqi are selected with 
an efficiency of more than 30% at our model points. Moreover, these efficiencies depend 
only weakly on the basic cuts on Ifi? and M e g. 

Evidence of squark pair production can be seen in the ratio of the events before and 
after the hemisphere cuts, because it is significantly different from that of gq production. 
At Point B, the ratio N(l + l + ; c 5 )/A(7 + Z + ; c ) is 0.171 with a statistical error of ~0.01. If 
qq production does not occur, N(l + l + ;co) becomes 1559, N(l + l + ;c^) becomes 80 and the 
ratio becomes 0.051 with a statistical error of ~0.006. Actually, there are various sources 
of systematic errors such as uncertainties on the PDF, NNLO corrections and so on. The 
total systematic error is 14% if the uncertainties from squark and gluino mass errors cancel 
in the ratio and we ignore errors on the branching ratios. However we think there is enough 
margin to identify qiqL production if careful analyses are done at LHC. 



IV. COMPARISON OF OTHER MODELS WITH THE MSSM 
A. The model with an extended gluino sector 

The model with an extended gluino sector has been discussed in section 2.2. In this 
model, a gluino acquires a Dirac mass term with an adjoint fermion a. As the majorana 
gluino mass parameter decreases from the MSSM value for the same gluino mass, the total 
SUSY production cross section decreases. In particlular, cr^g^) decreases more rapidly 
than <j{gq~L) when the majorana gluino mass parameter is reduced. Then SS2/ events from 
QlQl decrease more than those from gq^. Figure [TUh shows A(/ ± / ± ; Co) as a function of the 
majorana gluino mass. Here, we set the mass spectrum of this model as that of Point B 
and = rrig and m~ g2 = —3000 GeV, and branching ratios and efficiencies of cuts are the 
same as in the previous section. Moreover, we simplify our calculation by assuming that all 
events from qq production occur from q[q[ productions (see Appendix IA 3j) and that 



the non-g, non-g^ contribution does not depend on m g . We show N{l + l + ; Co) as a bold solid 
line, N(l~l~; Cq) as a bold dashed line, N(l + l + from qq; cq) as a thin solid line and N(l~l~ 
from qq; cq) thin dashed line. 
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Figure [TUb shows A^(/ ± / ± ;c5). SS2/ events from q^q^ are dominant in the total SS2Z 
events after the C5 cut, and that the dependency on the majorana gluino mass of the total 
number of SS2/ events is nearly the same as the dependence shown in Figure H On 

the other hand, the gq contribution is dominant in Figure [TUh . 

We show the ±lcr statistical error for the MSSM limits with a dark gray zone. The 
absolute numbers of SS2/ events after the cut cq, C5 depend on various parameters, acceptance 
and so on. We assume the total uncertainty is 17% as discussed in the previous section. 
This is also shown in Figure [10] with a light gray zone. 

The N(l + l + (l~l~); cq) shows more than 17% deviation from the MSSM in case of 
m g < 606 GeV (m g < 459 GeV). This means m D > 1262 GeV (m D > 1427 GeV). It 
is m g < 936 GeV (m g < 730 GeV) for N(l+l+ (/"/"); c 5 ). This means m D > 664 GeV 
(m D > 1089 GeV). 

As discussed in the previous section, the uncertainties on the leptonic branching ra- 
tios, PDF, QCD NNLO corrections and squark mass errors partly cancel by taking the 
ratio N(l + l + ;c 5 )/N(l + l + ]c ). Roughly speaking, iV(Z ± Z ± ; c 5 )/V(/ ± Z ± ; c ) depends linearly 
on a(q^q^)/a(gqL)- In this model, the ratio cr(q L q^)/a(gq^) decreases as the majorana 
gluino mass decreases from the MSSM value. The ratios N^^; cs)/iV(Z ± Z ± ; Cq) ctS ct func- 
tion of the majorana gluino mass are plotted for Points A~D in Figure [TTk~d. We show 
only the statistical error of ±ler in the MSSM limit for 3 x 10 5 events with a gray zone. 
For example, the ratio A^(/ ± / ± ; cs)/iV(/ ± / ± ; cq) has a statistically significant difference from 
Pont B for m g < 850 GeV at the la level. Because a(q~q~) <^ a(q + q + ) the sensitivity to 
is worse. 

B. The Littlest Higgs model with T-parity 

We now consider the case where the I$t signature arises from decays of the quark partner 
g_ in the LHT model. The g_ has a similar decay pattern to the q. Indeed, the leptonic 
branching ratio of qt is ~15%, which is almost the same as that of q + at Point C and D. 
The acceptance should be similar to that of the MSSM, because the acceptance depends on 
the decay kinematics, namely on the mass difference between q or g_ mass and the lightest 
R-odd or T-odd particle. 

While the collider signal is similar to that of qq production in the MSSM, the LHT 
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FIG. 10: Number of SS2Z events as a function of the gluino majorana mass at Point B a) after the 
cut Co, b) after the cut C5. Dark gray zones show la statistical errors for 3xl0 5 events. Light gray 
zones show 17% errors. 

model predicts different production cross sections from that of the MSSM. The cross section 
cr(qZq_ ) is 0.70 pb and a(qZqZ) is 0.15 pb for M q _ ~ 800 GeV. On the other hand, a(q^q^) is 
0.13 pb and a{q£q£) is 0.037pb at Point C (m 9 ~ ~ 800 GeV). For M q _ = 1000 GeV, <r{qtq±) 
is 0.22 pb and a(qZqZ) is 0.045 pb, while cr(q^q^) is 0.049 pb and cx(g^g^) is 0.012 pb at 
Point D (rriq ~ 1000 GeV). a(q^iq^) is about 4~5 times larger than c(q^q^) at each point. 
Note that cr(qiqL,) is always significantly smaller than o~(g_g_) of the LHT model no matter 
how heavy or light the g is. Moreover, there is no g production in the LHT model. If the 
excess of the production cross section is established and the existence of a light gluino is 
excluded, we can claim the LHT signature is observed. 

At Point C, N{l+l + ; c ) = 467, N(ri~; c ) = 233, N{l + l+; c 5 ) = 62 and N{l~r; c 5 ) = 17 
for an integrated luminosity of 86fb _1 . To study N(l + l + ) or N(l~l~), we need to simulate 
all LHT production processes, which is beyond the scope of this paper. We assume the 
number of signal events from production process is the number of signal events from 
qq production process scaled by the ratio of cross sections cr(q^q^)/a(qLqL), and the other 
production processes are ignored for simplicity. We expect N(l + l + ; Co) ~ 816, N(l~l~; cq) ~ 
162, iV(7 + Z+; c 5 ) ~ 302 and N(ri~; c 5 ) ~ 49 in the LHT model of M g _ ~ 800 GeV for the 
same integrated luminosity. 

At Point D, we find N(l + l+;c ) = 571, iV(rr;c ) = 263, N(l+l+;c 5 ) = 78 and 
N(l-l~;c 5 ) = 24 for an integrated luminosity of 307fb~ 1 . We expect N(l + l+;c ) ~ 962, 
N{1-1-]Cq) ~ 210, N(l + l + ;c 5 ) ~ 321 and iV(rr ; c 5 ) ~ 61 in the LHT model of 
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FIG. 11: N(c5)/N(cq) dependence on the majorana gluino mass m g at a) Point A, b) Point B, c) 
Point C and d) Point D. Gray zones show la statistical errors for 3xl0 5 events. 



M q _ ~ 1000 GeV for the same integrated luminosity. We can see that in the LHT models, 
N(l + l + ; Co) : N(l~l~; cq) would be about 4:1. This is different from the cases at the MSSM 
models Point C and D. 

If there is no particle production except g_ the ratio N(l + l + ; c 5 )/iV(/ + Z + ; c ) is expected 
to be around 0.37 in the LHT model of M q _ ~ 800 GeV, which is larger than the value 
at Point C (0.133). The ratio N(1 + 1 + \c§)/N(1 + 1 + \cq) is expected to be around 0.33 for 
M q = 1000 GeV, and the value is larger than that at Point D (0.137). 

Finally we comment on the case where q and g are highly degenerate so that we cannot 
detect the jets from g — > qq decay by hemisphere analysis. Note that cr(gq) ^ a (QQ)i 
therefore the rate of the SS2Z events could be as large as the LHT prediction. However, even 
in this case, N(l + l + ; cq) : N(l~l~; cq) in the MSSM is ~ 2:1 and cannot be 4:1 because there 
are gq and gg contributions. 
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V. CONCLUSION 



Information on the fundamental Lagrangian of SUSY model can be extracted from each 
SUSY production process at LHC. For example, q^L production cannot occur without a 
majorana gluino mass, because a chirality flip is required. Thus the majorana nature of 
gluino mass can be extracted from this process. 

This process can be investigated using the SS2/ events because B1Z(qf — > + X) ^> 
B1Z(qf — > l T + X). At LHC, however, mixed production with other sparticles g, q, etc. 
make it difficult to interpret the signal. In particular, gq production also contributes to the 
SS2/ channel. 

In this paper, we have discussed a systematic method to separate the production modes. 
When we measure a(qiqL — > l^l^ + X) in the MSSM, we suffer from a problematic back- 
ground from gqi production. We have proposed a new method based on a hemisphere 
analysis as a solution to this problem. In the hemisphere analysis, we assign high p T objects 
into two hemispheres, where each hemisphere contains high pt objects from the same parent 
particle with high probability. Then we require that there is only one jet with p T > 50 GeV 
in a hemisphere. For the sample after some basic cuts, 30~40% for q L q L and 1~3% for gq L 
remain after the hemisphere cut. Therefore we can obtain SS2/ events with enhanced qiqL 
contribution, which may be used to estimate oi^q^qi). 

We have also discussed two models which have similar collider signals to the MSSM but 
whose relevant production cross sections are dramatically different. 

One of the models is the MSSM with an extended gluino sector, where the gluino can 
have a Dirac mass with an adjoint fermion a. oi^q^qj)) is sensitive to the fraction of majorana 
mass terms in the gluino mass. In the case where the gluino is pure Dirac, q^L production 
cross section becomes zero. 

We have applied our analysis to the model points with an extended gluino sector which 
have the same mass spectra of some MSSM model points except for an additional heavy 
adjoint particle. We estimate the number of SS2/ events as a function of the majorana gluino 
mass, and estimate the sensitivity to the Dirac gluino mass. We take only the statistical 
error into account and assume that the masses of the squarks and the branching ratios are 
known. We find that the ratio a(q^q^) /ar(gq^) that can be estimated from the acceptance 
under the hemisphere cut is useful because this quantity should be less sensitive to the error 
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on the parameters and the acceptance and uncertainties of the PDF and QCD corrections. 

We have also considered the LHT model. In this model, a set of T-odd partners is 
introduced to the SM matter particles and the EW gauge bosons. The quark partner 
production cross section cr(q_q_) is 4~5 times as large as the production cross section 
^{QlQl) in the MSSM. The g_ and q have similar decay patterns and branching ratios. 
Thus, the number of SS2/ events from g_g_ productions in the LHT model is expected 4~5 
times larger if patterns of the mass spectra are the same. Moreover, the LHT model has no 
process corresponding to the process gq nor gg production in the MSSM. To exclude gluino 
productions is important to identify the LHT model. This can be done by investigating 
N(l + l+)/N(l-l-). 

It is generally important to measure the production cross sections of sparticles separately 
to verify the MSSM and distinguish various models. In this paper, we develop a method to 
identify gluino and squark production separately for SS2/ channel. The method is based on 
the cuts on the kinematical configulation of the jets and can be applied to the other models. 
More development is needed for the model independent study of physics beyond the SM. 
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APPENDIX A: APPENDIX 



1. Mass spectra at our model points 

We show mass spectra for the selected model points which have been analyzed. They are 
calculated by IS A JET 7.72. They are all mSUGRA mass spectra except that gluino masses 
of Points C and D are 300 GeV larger than mSUGRA predictions. 
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mass parameter 
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9^fi 3fi 


394 47 


437 1 3 
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e R 
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3Q3 OQ 


43^ 


1 43 00 


Ve 


243.67 


314.44 


429.70 


496.84 


186.00 


T\ 


160.79 


194.28 


387.18 


429.40 


133.39 


Tl 


258.50 


325.66 


437.28 


503.05 


206.02 


V T 


242.89 


313.50 


428.00 


495.00 


185.06 


X? 


132.74 


179.11 


133.95 


180.54 


96.05 




253.87 


345.39 


256.39 


348.15 


176.80 


xl 


-448.66 


-575.60 


-451.46 


-577.05 


-358.82 


xl 


467.28 


591.37 


470.52 


592.66 


377.84 


xt 


254.13 


345.95 


256.71 


348.72 


176.37 


xt 


466.48 


590.74 


469.74 


592.56 


378.26 



TABLE X: Mass spectra of sparticles for the selected model points 
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2. Branching ratios at our model points 

We also show the branching ratios of the selected points. They are also calculated by 
ISAJET. Squarks mainly decay into charginos X s1 and the second lightest neutralino x° & t 
each point. Note that, at Points C and D, x° does not decay into sleptons. On the other 
hand, at Points A, B and SPSla, x° decays into or Xi T+T ~ ■ 
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1.4 


1.4 


Xi~ 




21 


7.7 


100 


100 


1.1 


in 




38 


36 


33 


38 


43 






24 


54 








0.4 






24 


28 


31 


37 







— ► X^T~i> T 


54 


38 








98 






10 


12 


12 





14 


X°2 




45 


55 


100 


100 


0.8 






24 


21 


20 


22 


36 




/v, 1 


7.6 


18 








13 




xr^ 


15 


14 


4.9 


4.4 


21 




X?T+T- 
/v. 1 


46 


20 








87 




— > X2 t 


34 


35 


39 


41 





xi 




14 


13 


14 


14 


13 




-► 


13 


13 


13 





35 




^x° 2 x 


24 


26 


25 


27 


23 






19 


21 


23 


31 


15 




-^xfw T 

/Y 1 


29 


29 


30 


30 


60 






9.6 


8.3 


3 


2.5 


17 


xl 




21 


19 


12 


13 


8.5 




-X> 


3.7 


3.5 


7.9 


9.4 


5.4 






18 


22 


22 


25 


15 






5.6 


5.2 


9.6 


11 


7.4 




^xfWT 


26 


26 


32 


30 


52 



TABLE XI: Branching ratios of sparticles for our model points: here, X means some SM particles. 
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3. Dominant processes contributing to SS2Z events 

Here we show the numbers of SS2Z events from each production process at Point B. Note 
that SS2/ (l^l^) signals come not only from qfqf production but also from q^ql (mainly 
ul(Il) production. For l + l + events, about 6% of total l + l + events are from qq production. 
The / + / + events from qiql production are 10% of total l + l + events from qq production. 
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Point B 


all 


CO 


Cl c 2 


c 3 c 4 


C5 ratio 


1+1+ 


QQ 


1479 


1014 1001 894 363 883 361 0.356 






1161 


788 


778 701 291 692 289 0.367 




ulcl 


95 


68 


68 58 


22 58 


22 0.324 






84 


58 


56 51 


17 50 


17 0.293 




u L d*r 


90 


62 


61 54 


19 53 


19 0.306 




U L S*r 
Li 


48 


37 


37 29 


13 29 


13 0.351 


1+1+ 


gq 


1765 


1098 


433 613 


52 279 


36 0.033 






1581 


974 


388 538 


41 252 


30 0.031 




gcL 


22 


18 


8 8 


2 3 


1 0.056 






32 


22 


8 11 


2 5 


1 0.045 




~gd* L 


78 


53 


17 38 


6 12 


3 0.057 




~ ~* 
9 S L 


39 


22 


9 13 


6 






Point B 


all 


CO 


Cl c 2 


c 3 c 4 


C5 ratio 


rr 


QQ 


519 


333 


323 293 


123 285 122 0.366 






280 


188 


181 164 


74 159 


73 0.388 




dhSL 


87 


62 


61 56 


31 55 


31 0.5 




ULdh 


106 


67 


65 57 


16 55 


16 0.239 




d L u* L 


20 


7 


7 7 


7 







d L c* L 


10 


7 


7 7 


1 7 


1 0.143 


rr 


M 


881 


553 


241 310 


23 172 


15 0.027 




g d L 


615 


396 


184 232 


21 132 


14 0.035 




g s L 


30 


18 


8 9 


5 







g u L 


135 


75 


23 40 


2 17 


1 0.013 




g u l 


67 


43 


19 21 


13 







— 

g c l 


20 


12 


4 6 


3 






TABLE XII: Contributions from each production process to SS2/ events at Point B. 
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